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Abstract 

The short distance behaviour of massive integrable quantum field theories is ana- 
lyzed in terms of the form factor approach. We show that the on-shell dynamics is 
compatible with different definitions of the stress-energy tensor T flll {x) of the the- 
ory. In terms of form factors, this is equivalent to having a possible non-zero matrix 
element F\ of the trace of T^ u on one-particle state. Each choice of F± induces a 
different scaling behaviour of the massive theory in the ultraviolet limit. 



1 Introduction 



One of the most relevant steps towards the physical interpretation of a given field dynam- 
ics consists in the identification of the stress tensor T^ix). It gives the local distribution 
of energy and momentum, and rules the response of the system under local scale trans- 
formations. In this paper we consider quantum field theories involving one scalar field 
(f(x) defined in (1 + 1) flat space-time where there is, however, a natural ambiguity in 
the definition of T liv {x). To see this, the simplest way is to initially consider the theory 
defined on a curved two-dimensional manifold with metric tensor g^ u and scalar curvature 
R 

A = J d^^Tg^d^d^ - v(<p) + aR^ , (i.i) 

In the above action, a is a free parameter that in the limit of a flat manifold g^ v — > 
r/^, labels the one-dimensional family of stress-energy tensors associated to the on-shell 
dynamics 

T^{x) = f^{x) + a{d^d u - g^n)(p(x) , (1.2) 



where 



T^(x) = 2tt 



;i.3) 



Although in the flat Minkowski space the equations of motion for the field ip do not 
depend on the specific definition of T^ u (x), the ultraviolet behaviour of the theory may be 
extremely sensitive to any modification of this operator. A well-known example is provided 
by the conformal invariant QFT (l|, [|: in a Coulomb Gas language, an extra derivative 
term in the definition of the stress-energy tensor results in a non-zero background charge 
that, in turn, induces a non-trivial scaling behaviour of the operators of the theory 0, |[]. 
It is therefore an interesting question to see what are the consequences of a redefinition 
of T^ v (x) also in the case of quantum field theories which do not have scaling invariant 
properties. A special class of these theories are the massive integrable models where 
the off-shell dynamics may be completely characterized in terms of the on-shell scattering 
amplitudes of the massive excitations. In fact, we can take advantage of their integrability, 
and compute exactly the matrix elements of local operators Ok(x) on the asymptotic 
states, by means of the Form Factor Bootstrap Approach [5-13]. All correlation functions 
are then reconstructed in terms of their spectral representation. For large values of the 
relative distances, these correlation functions will have an exponential decay ruled by 
the lowest massive state appearing in their spectral decomposition. Their short distance 
behaviour, on the contrary, present power law singularities. In order to identify the 
conformal dimensions of the operators, we need to analyze the short distance singularity 
of the correlation functions < T zz (z,z) 0^(0) >. For those operators that correspond to 
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primary fields in the conformal limit (mR) — * (R 2 = zz), we have 

<T zz (z,z)O k (0)>~ . (1.4) 

The conformal dimensions A(O k ) will depend in general on the definition of T liV {x) and 
therefore different definitions of this operator may induce different scaling behaviours of 
the theory in the deep ultraviolet region. Analogously, the central charge of the ultraviolet 
theory may be extracted from the short- distance singularity of the correlator 

<T zz (z,z)T zz (0)>~ JL . (1.5) 

The plan of the paper is as follows. In section 2 we discuss the general properties of 
the stress-energy tensor. Using entirely the formalism of the form factor approach, we 
derive the conditions on the matrix elements of T^x) which guarantee its conservation 
and locality. From this analysis, there is a one- dimensional space of possible stress-energy 
tensors T^ u (x) compatible with the on-shell dynamics. In section 3 and 4 we then discuss 
the ultraviolet behaviours of the simplest Afflne Toda Field Theories p3] , i.e. the Sinh- 
Gordon and the Bullough-Dodd models based respectively on the simply-laced algebra 
Ai and on the non-simply laced algebra A% ■ In both theories, different ultraviolet 
limits may be reached by varying the definition of the stress-energy tensor, with the 
relevant conformal data simply obtained in terms of the form factors. Our conclusions 
are summarized in section 5. 



2 Form Factors of the Stress-Energy Tensor 

Let us consider an integrable two-dimensional massive QFT characterized by its elastic 



factorizable S-matrix [ T%, [H| [L6| . As it is well known, the form factor approach is quite 
effective to characterize the operators in such a theory in terms of their matrix elements 
on the set of asymptotic states ]5|, [|. Adopting the standard parametrization of the 
momenta in terms of the rapidity variable (3, i.e. = (m cosh/5, m sinh/?), any local 
operator Ok(x) will be uniquely identified^] by the set of its Form Factors (FF) 

F n fc (/3 ...,/?„) =<O|0 fc (O)|/3i,...,&>to. (2.1) 

In the above definition, the set is ordered as ft > j3 2 ■ ■ ■ > Pn and our normalization 
is given by 

n 

in <#,..., (3' m | ft, . . . , (1 n > in = 8 m>n J] 27r5(# - A) • (2.2) 

8=1 



1 This justifies the interchangeable use we make of the words "operator" or "form factors" in the rest 
of the paper. 
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Once the FF of the local operators Ok(x) are known, their two-point (and higher) corre- 
lation functions can be reconstructed through the unitary sum0 



< O k (r) O k (0) > = £ f d/h U0 -^ n I F*(A ...p n ) | 2 exp (-mr £ cosh t ) (2.3) 



Let us discuss then the properties of the FF, as dictated by relativistic invariance and 
general requirements of QFT H [|. For an operator O k (x) of spin s, relativistic invariance 
implies 

F*(/3 1 +A,ft + A,...,/3 ft + A) = e sA F*(p u (3 2 , . . . , p n ) . (2.4) 

The FF of a given theory are solutions of a set of functional and recursive equations. 
The functional equations arise from the monodromy properties of the functions F% which 
are ruled by the S-matrixP] 

F*(/V--,A,/W--,/y = F^p 1 ,...,p i+1 ,p i ,...,p n )S(p i -p i+1 ) , (2.5) 

n 

F^p 1 + 2m,...,p n _ 1 ,p n ) = F n fc (/5 2 ,...,/3 n ,/3 1 ) = n5(A-/3i)F n fc (/5 1 ,...,/3 n ) . 

The recursive equations are obtained, on the contrary, by looking at the pole structure of 
the matrix elements F%. The first kind of poles are kinematical poles located at (3^ = in. 
The corresponding residues give rise to a recursive equation between the n-particle and 
the (n + 2)-particle FF 

-^lim(/3-/?)F n fc +2 (/3+^7r,/3,/3 1 ,/? 2 ,...,/3 n )= fl - f[ S{(5 - ^)) F^ u ...,f3 n ) . (2.6) 
0-0 V i=i / 

If bound states are present in the spectrum, there is another set of recursive equations 
obtained by looking at their poles in the matrix elements. Let /3y = i u be the location of 
the pole in the two-particle scattering amplitude corresponding to the bound state. Then 
the corresponding residue in the FF is given by 

-zlimeF n+1 (/3 + z«-^/3+^/3 1 ,...,/3„_ 1 ) = FF n ((3, ft, . . . , (3 n _ x ) , (2.7) 

where T is the on-shell three-particle vertex. This equation establishes a recursive struc- 
ture between the (n + 1)- and n-particle form factors. 
The two chains of recursive equations, 

• • • — > F n+ 4 — > F n+2 — > F n — > i^„_2 — * ■ ■ ■ (2 8) 

. . . — > - -Pri+3 — * F n+ 2 — > F n+ i — > . . . 



2 This expression holds for scalar operators. The generalization of this equation for operators of spin 
s is easily done by using eq. (2.4). 

3 For simplicity we consider theory with only one massive self-conjugate particle. 
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(and the consistency conditions associated to them) are quite effective for the explicit 
determination of the FF of a given theory. 

The above discussion holds for any FF of a local operator, in particular for those of the 
stress-energy tensor T^ v (x). However, due to the special role of this operator, its FF has 
some distinguishing properties. Since T^ix) is conserved, it may be expressed in terms 
of an auxiliary scalar field A(x) as 

T^(x) = (d^d u - g^U) A(x) . (2.9) 

In light-cone coordinates x ± = x° ± x 1 , its components are given by 

T ++ = d\A ,T__ = d 2 _A , (2.10) 

Q = T^ = -DA = -4d + d_A . (2.11) 
Introducing the variables Xj = and the elementary symmetric polynomials in n- 

(n) 

variables a\ defined by the generating function 

n n 

no* + *j) = £^, (n w ••,**) . ( 2 - 12 ) 



j=l i=0 

it is easy to see that 



1 , la^' 2 



F n ^(/3 1; ...,/3 n ) = - i mM^l F^(p 1 ,...,p n ) , 

Fj--(A,-..,/3„) = ~\m 2 (^ ) ) 2 F n A (/3 1 ,...,/? n ) , (2.13) 

(«•)(«) 

if (/?!,...,/?„) = m*^^±F A ((3 u ...,(3 n ) . 



Solving for F A , we have 



1 a [n) 



Fn ++ (Pi, ■ ■ ■ , Pn) = ~ ~ -J^TJn) F n(Pl, •••,#») 

i W (») 

Fj-(A,...,/? n ) = - 7 ^^F n e (/5 1 ,. ..,/?„) • (2-14) 



Hence the complete knowledge of is encoded into the form factors of the trace 9. 
As any spinless operator, its form factors F®(/3i, . . . ,/3 n ) depend only on the difference 

of the rapidities = $ — Moreover, since the FF of T and T ++ have the same 

singularity structure of the FF of 6, F®(/3i, . . . , /3 n ) (for n > 2) has to be proportional to 
the combination of symmetric polynomials Oi<Jn-i which corresponds to the invariant 
total energy-momentum. 
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Additional constraints on are obtained from the knowledge of their asymptotic 
behaviour in each variable ft. This behaviour generally depends on the particular model 
under consideration. For the lagrangian quantum field theories discussed in this paper, 
we have 

F n e (ft + A,ft,...,ft) A -=5? (i) , (2.15) 

i.e. they become constant for large values of the individual momenta. This condition 
can be easily checked by analyzing the asymptotic behaviour of the Feynman diagrams 
entering the perturbative definition of these matrix elements [[], 

Concerning their normalization, the recursive structure of the space of FF reduces 
the problem of finding the normalization of the matrix element of Q(x) for the initial 
conditions of the double chain (|2.8|) , i.e. the two-particle FF -F 2 (ft 2 ) and the one-particle 
FF F*(J3). 

The normalization of the two-particle FF F®(f3i 2 ) can be fixed by making use of the 
definition of the energy operator || 

1 r+oo 

E = — dx 1 T 00 (x) . (2.16) 

2% J-oo 

In fact, computing the matrix element of both terms of this equation on the asymptotic 
states < ft| and \fi >, for the left hand side we have 

< f3'\E \(3 >= 2itm cosh f3 5(p' - (3) . (2.17) 

On the other hand, taking into account that T 00 = d\A and using the relation 

< ft| 0(x)\f3>= e ^^- pM (^F 2 °(/3, ft - in) (2.18) 

valid for any hermitian operator O, we obtain 

F? A (Pi, ft) = - m 2 (sinh ft + sinh/3 2 ) 2 F 2 A (P 12 ) . (2.19) 

Then, from eqs. ( 2.13 ) and ( 2.16|) , the normalization of F 2 @ is given by 

F 2 e (in) = 2nm 2 . (2.20) 

However, no special constraint exists for the matrix element of Q(x) on the one-particle 

state 

F® =< | 6(0) | > , (2.21) 

which is then a free parameter of the theory. Notice that from Lorentz invariance, it does 
not depend on the rapidity variable ft. 

Since higher FF of G are obtained as solutions of the recursive equations ( |2.6|) and 
( |2.7| ) (with initial condition given by the one-particle and two-particle matrix elements), 
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the arbitrariness of F® propagates in the recursive structure (2~E) of the FF and therefore 
gives rise to a one-parameter family of possible stress-energy tensor T^ u for a given theory. 

A simple example of the above discussion is provided by the free massive bosonic 
theory, with equation of motion 

(□ + m 2 )y? = . (2.22) 

The 5"-matrix in this case is simply S = 1 and therefore the FF have trivial monodromy 
properties and simple analytic structure. Among the local operators of the theory, the 
elementary field ip(x) is identified by the set of FF 

W^.-A) =<0b(0)|ft, &,...,&>= ±=5 x , n , (2.23) 

and its two-point euclidean correlator reduces to a Bessel function 

^Tq Til J-oo Z7T J-oo lH 

= -K (mR) . (2.24) 

7T 

The absence of interaction implies that the composite operators ip k /k\ are simply defined 
by the following FF 

< 0|^|/3i, /V-., P n >= (^j 5 n , k . (2.25) 

The equation of motion is compatible with a class of stress-energy tensor labelled by the 
free parameter Q appearing in the definition of Q(x) 

Q(x) = 2ir (mV + ^ D ^) • ( 2 - 26 ) 

In terms of FF we have 



e 

Ff = -V27rm 2 Q , (2.27) 



F- = 



F 2 e = 2vrm 2 , 
Ff = ,k>2 . 



The meaning of Q becomes clear once we analyze the ultraviolet limit of this massive 
theory. The central charge of the underlying CFT which rules the ultraviolet properties 



of the model may be computed by using the c-theorem sum rule |L8|, nj 

3 r°° 



poo 

/ dRR 3 < Q(R)e(0) > E ■ (2.28) 
Jo 
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The euclidean correlator is given by 

<Q(R)Q(0)> E = fy^^..J^f^\F k (pi,...,p k )\\-- R ^^ 

= m 4 (2 (K (mR)) 2 + 2Q 2 K (mR)) , (2.29) 
and the result of the integral ( 2.28Q is simply 

c = 1 + 12 Q 2 . (2.30) 

Hence the one-particle FF of G is related to the background charge of the CFT reached 
in the ultraviolet limit. We could have obtained the same result by directly analyzing the 
ultraviolet limit of the holomorphic component of the stress-energy tensor. Indeed 

^2 



<T zz (z,z)T zz (0,0)> E = (-J <T zz (R,R)T zz (0,0)> E 

= m 4 (2 (K 2 (mR)) 2 + 2Q 2 K 4 (mR)) , (2.31) 

(R 2 = zz and z = x° + ix 1 , z = x° — ix 1 ), and in the limit (mR) — > we have 

<T zz (z,z)T z M0) >= = 1 \l 2 4 Q2 ■ (2.32) 

To complete our discussion on the free theory, let us compute the conformal dimensions 
A (a) characterizing the scaling properties of the exponential operators V a = e aLp in the 
ultraviolet limit. This will be identified as the coefficient of the most singular term in the 
ultraviolet limit of the correlator < T zz (z, z) V a (0) >. Using the FF of V a (0) given by 



/ 



n 



<0\V a (0)\p 1 ,...,p n >= \-j=\ , (2.33) 



we have 



<T zz (z,z)V a (0)> E = (-) <T zz (R,R)V a (0)> E 



and for (mi?) — > 



m 2 (-^-( Kl (mR)) 2 + a-Q=K 2 (mR) 



< T zz (z, z) V a (0) > E ~ \ % (-^ + ^| ) , (2.34) 



i.e. 

A(a) = — + ^=. (2.35) 

Stt V47T 



Due to the background charge, they differ from the gaussian value A(a) = —a 2 /87i. 
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3 The Sinh-Gordon Model 



The free theory provides an easy example of non-trivial ultraviolet behaviours induced by 
a non-zero value of Ff. It is interesting to see if similar occurrences are also present for 
interacting theories. In this section we will analyze the Sinh-Gordon model and in the 
next one the Bullough-Dodd model. 



3.1 Basic Properties 

The Sinh-Gordon theory is a classical integrable model defined by the equation of motion 



2 

U(p = (e 9tp - e~ g<p ) . (3.1) 

The theory is invariant under a Z 2 -symmetry ip — > —(p. The integrability of the model 
also persists at the quantum level and the exact two-body elastic S'-matrix involving the 
asymptotic particles created by the operator ip is given by [T7| 

o (B m _ tanhj(^-if) 

S( ^ B) ~ tanh±(/? + ? f) ■ (3 ' 2) 
The coupling constant dependence of the model is encoded into the function 

Q 2 1 

Big) = — — . (3.3) 

yyj 47r(l+^ 2 /87r) v ; 

Since there is no pole on the physical sheet < Im (3 < n for real values of the coupling 
constant g, the Sinh-Gordon model presents no bound states. Notice that the S'-matrix is 
invariant under the duality transformation B — > 2 — B, which establishes a mapping of 
the theory between the weak coupling and strong coupling regimes, i.e. g — > 8ir/g. This 
symmetry will be respected by all FF of manifestly self-dual operators. 



3.2 Space of the Form Factors 



The form factors of the Sinh-Gordon model have been investigated in [|| |10|. In this 
subsection we briefly recall the basic results obtained in refs. P, [T(| which are relevant for 
our subsequent considerations. 

In order to compute the FF of this theory, the first step is to take into account their 
monodromy properties dictated by the S'-matrix. To this aim, let F^* n (/3,B) be the 
solution of the functional equations 

5) = B) Si/3, B) , 

F^n-(3,B) = F^n + ^B) , 
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which has no poles and zeros in the physical sheet < Im (3 < it and with an asymptotic 
behaviour given by 

(3.5) 



jiniC(AB) = l 



Its explicit expression reads 



F^{p,B) = Af(B)exp 



sinh (if ) S inh(f(l-f)) sinhf 



o x 



JV(B) 



exp 



sinh x 

rf x sinh(f ) sinh (§(!-§)) sinhf 



sm 



'x£ 

2?r 



X 



sinh x 



(3.6) 



where (3 = in — j3. Equivalently 



F*tw,B) = n 



fc=0 



r(^ + l + S)r(^ + l-f + I)r(^ + i + f + f 



(3.7) 



Since S(0,B) = —1 (for B 7^ and 2), F^ n ([3,B) vanishes at the two-particle threshold 



value (3 = 0. As discussed in |TT|, this generally induces a suppression of all higher 
thresholds in the spectral representation of the correlation functions and gives rise to very 
fast convergent series. 

In terms of F^ n (f3,B), a convenient parameterization of the n-particle FF of the 
Sinh-Gordon model is given by 



F n (/3i,...,A0 = H n Q n (x u ...,x n ) n ?^''? 

i<j \ X i + X j) 



(3.8) 



Here Xi = e^' and if n are normalization constants, which can be conveniently chosen as 

H 2n+1 = H lf i 2n (B) , H 2n = H 2 fi 2n ~ 2 (B) , (3.9) 



with 



/ 4sin(7ri?/2) \ 



(3.10) 



The functions Q n (xi, . . . , x n ) are symmetric polynomials in the variables x%, . . . , x n , solu- 
tions of the recursion equations 



( — ) Qn+2^~ x i x i x li ■ ■ ■ 1 x n) — xD n (x, X\, X 2 , ■ ■ ■ , X n ) Q n [ x lj x 2i ■ ■ • > x n) > (3-H) 



with 



Dn\X) X\j . . . , x n j 



n 

E 

k=l m=l,odd 



E 



2(n-fc)+m (n) (n) 



a k a k-m\ 



-1) 



fc+1 



(3.12) 
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and 



sin(mrTr) 



N = — ( i; . (3.13) 

sm(7r -j) 



For FF of spinless operators, their total degree is equal to n(n — l)/2 whereas their partial 
degree in each variable Xi is fixed by the asymptotic behaviour of the operator Ok which 
is under investigation. 



As shown in the problem to classify all possible scalar operators of the Sinh- 



Gordon model reduces to find the most general class of solutions of eq. ( |3.11| ). Since this 
is a homogeneous equation, its solutions span a linear space whose basis may be written 
in terms of the so-called elementary solutions given by 

Qn(k) = detMy(fc) , (3.14) 

where Mij(k) is an (n — 1) x (n — 1) matrix with entries 

Mij(k) = a 2i -j + k] . (3.15) 
These polynomials depend on an arbitrary integer k and satisfy 

Q n (k) = (-l) n+1 Qn(-k) . (3.16) 

Therefore the structure of the FF of the SGM consists in a sequence of finite linear spaces 
whose dimensions grow linearly as n increasing the number 2n — 1 or In of external 
particles. The reason is that, at each level of the recursive process, the space of the 
FF is enlarged by including the kernel solutions of the recursive equation (|3.11|) , i.e. 
Q n (—x, x, x\, . . . , x n -2) = 0. With the constraint that the total order of the polynomials 
is ra( - n ~ 1 ^ , the kernel is unique and given by E n (xi, . . . ,x n ) = det oii-y These solutions 
then gives rise to the half-infinite chains under the recursive pinching x\ = —x<i = x 



Qn+4 ~ * Qn+2 ~ > Qn ~ > Qn-2 ~ * • • • ~~ * Q3 ~ * 1 

_ 2 — > ■ ■ ■ — > ^2 



y y Qn ^ ~~ > Qn- 



n (n-2) _ n {n-2) _ ^(n-2) 

Vn+4 ' Vn+2 ' V n 

n( n ) _ /)(") _ v 

Vn+4 — ' Vn+2 " ' 

Q 



(3.17) 



n+4 ~~ ' Vn+2 
(n+2) _ v 
^4 ' ^n+2 



The explicit expressions of such solutions can be found by determining the linear combi- 
nation of Q n (k) which reduces to S n at the level n. 
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3.3 Cluster Operators and Fundamental Exponentials 

The fundamental exponential operators &±(x) = e ±gip ^ define the Sinh-Gordon model 
and in general appear in the expression of the stress-energy tensor. In order to calculate 
their matrix elements, let us consider initially those FF which satisfy the requirements 

• To be asymptotically constant for /3j — > oo, i.e. 

F n ((3x + A, (3 2 , ...,P n ) A -=^ o(l) . 

• To be proportional to the combination?] o\o n -\ (for n > 2). 

• To be the solution of the cluster equations 

Jim F k+l (ft + A, . . . , fa + A, Pk+t, pk+i) = F k (/3i,..., Pk) F (p k+l , p k+l ) 

with initial condition Fq = 1. 

There are two classes of FF which fulfill the three above conditions. Their expressions 
are given by 

= Ht\B)Q n (l) n F f n( ^'f , (3.18) 

i<j \ x i i x j) 

where 

H£\B) = MB))" , HP = (-1T(KB)T ■ (3-19) 

The corresponding operators, which are self-dual by construction, will be called cluster 
operators and denoted as V±(x,B). We conjecture that the fundamental exponentials of 
the Sinh-Gordon model may be written asQ 

$ + (x,B) = B{l-B) V + {x, B) + 8{B-1) V-{x,B) , 
$_(x, J B) = e(l-B)V-{x,B) + B(B- l)V+(x,B) . 

Postponing a non-trivial check of this position until when we will study the UV-behaviour 
of the model, let us in the meantime discuss the properties of the operators $±(x, B) so 
defined. 

First of all, they satisfy the cluster property by construction, in agreement with the 
perturbative analysis for the matrix elements of the operators e ±g '~ p ( x \ Secondly, the 
FF of <&±(x,B) are not individually invariant under the duality transformation but each 
operator is mapped onto the other under the mapping B — > 2 — B, i.e. 

$±{x,B) = $ T {x,2-B) . (3.21) 



4 As discussed in sect. 2, this factorization property is shared by the general FF of 9. 
5 For the value of the step function at the origin we use 8(0) = 1/2. 
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Therefore they form a bidimensional representation of the duality symmetry. However, 
this mapping becomes degenerate at the self-dual point B = 1 where an identification 
occurs between the two exponential operators <&±(x,B). Namely, at B = 1 the matrix 
elements of the two fundamental exponentials become indistinct and denoting by $(x) 
the resulting operator, its FF are given by 

F$(/3i = f (M!)) n <3n(l) nr<i^fn(/3y)/(^ + ^) neven, 

I n odd . 

The identification of Q± (x) at the self-dual point has the additional consequence that the 
resulting field is an even operator under the Z2 parity of the Sinh-Gordon model, as 
it is evident from the vanishing of its matrix elements on all 2n + 1 particle states. 

Using the FF of the fundamental exponentials and those of the elementary field given 

by 

F*{p u ...,p n ) = H%Q n (0) n mi t J) (3-23) 

HL + i = ^= (KW , Ht n = , 
the quantum version of the equation of motion may be written as 

n<p(x) = (9(1 -B)- 9(B - 1)) ( e -«*0 - e^)) 

9 

771 

^ ^Mb) {v - {x - b) - vm ■ (3 ' 24) 

This equation has to be understood as an identity satisfied by the FF of the operators 
appearing on the left and right sides of this relation. 

3.4 Class of Stress-Energy Tensors 



The quantum equation of motion ( 3.24 ) is compatible with a one- dimensional space of 
stress-energy tensors given by 



Q(x) = Fg{B) {a$+(x,B) + (1 - a)$_(x,B)) . (3.25) 
The normalization constant F®(B) may be fixed by means of the Thermodynamical Bethe 



Ansatz 20 



^ - i^Sm ■ < 3 - 26 > 

The variable a, on the contrary, is a free parameter. Varying its value, we can weight 
differently the two fundamental exponentials in the trace and, consequently, we can in- 
terpolate between different scaling regimes of the Sinh-Gordon model in its ultraviolet 
limit. 
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3.4.1 The case a = 1 



The trace of the stress-energy tensor is given in this case by 

e(a;) = F°(B)$+(x,B) . 



(3.27) 



With such definition of 0, we expect that the massive theory will flow in the ultraviolet 
regime to a CFT defined by the bare action 

,2 



5_ 



d x 



2 lW 2g 2 



(3.28) 



In order to support this conclusion, let us compare the central charge associated to the 
CFT ( |3.28| ) with the central charge obtained, on the contrary, in terms of the FF by using 
the c-theorem sum rule Q2.28 ). 

Assuming eq. ( |3.28| ) as definition of the ultraviolet theory, the corresponding central 



charge is given by (see, for instance [21 



c(g) = 1 + 12 Q 2 .^) 



where 



Q~(g) = - 

Using eq. (|3.3|) , it may be written as 

c(B) = 1 + 6 



An 



B 



+ 



9 



2V47T 



B 



B 



B 



(3.29) 



(3.30) 



(3.31) 



Notice that this is a self-dual function of the coupling constant, i.e. invariant under 
B ->2- B. 

On the other hand, we may compute the central charge associated to the ultraviolet 
limit of the massive theory in terms of the second moment of the two-point function of 
the trace G(x) [|TJ, 

3 



1 f°° 

c(B) = - / dRR 3 < 0(i?)0(O) > E 

2 jo 



(3.32) 



According to eq. (|3.27 ), the two-point function of the trace 0(x) has to be computed in 
terms of the FF of the operator $ + {x, B) defined in eq. ( |3.20| ). The data reported in Table 
1 and plotted in Fig. 1 show that the first two FF of <&+(x) are sufficient to saturate the 
sum-rule ( [3.32 ) and to reproduce with high percentage of precision the expression ( 3.31 



An additional check that the ultraviolet limit induced by this choice of is ruled by 
CFT (|3.28|) , is given by the computation of the conformal dimensions of the fundamental 
exponentials. This can be done in two different ways, using directly CFT method or 
analyzing the ultraviolet behaviour of massive correlators. 
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For the CFT defined by eq. ( |3.28| ), the conformal dimensions of the primary fields 
corresponding to the operators e a(p are given by 

A -(") = -f- + — 7j=~ > (3-33) 

07T V47T 



and therefore, for the fundamental exponentials we have 

A -'*-) = 1 ■ . (3.34) 

On the other hand, we may compute the conformal dimensions of the fundamental expo- 
nentials by investigating the UV-limit of the correlators 



n=0 

x 



Z rZn n! J-oo 27T J-oo 2tt v 

F*±(/51,...,/? n )e- miJ ^ coshft . (3.35) 



At order 0(g 4 ), it is sufficient to truncate the series to the first two terms and also use 
the perturbative expansion 

,2 



M{B) = 1 - ^ + 0(g*) . (3.36) 



Since for small values of g 



(n) 

F^( Xl ,...,x n ) = -F°(B) Fp( Xl ,...,x n ) , (3.37) 



and 



<T„(z,z)$±(p)> E ^ -\-F®(B)F^F^^-K 2 (mR) (3.38) 

2 I 47T 



2?r 2 

in the limit mR — > 0, we have 



d/5 K 2 (2mR cosh /3) F 2 y+ (2/3)F 2 y± (2/?) 



A _($+) 

< T ZJS (z, Z) $±(0) > E ~ ^ > (3-39) 

z z 



with 

A_($_) = l + 0(<? 4 ) 



(3.40) 



in agreement with eq. Q3.34 ) . 
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3.4.2 The case a = 



The trace of stress-energy tensor is given in this case by 

Q(x) = F®(B)$4x,B) . (3.41) 

and we expect that the massive model will flow in the ultraviolet limit to a CFT defined 
by the bare action 



S+ = d x 



?- 


2 

m 

e 




2<? 2 


by 







9 \ 



(3.42) 



Q+(9) = (— + 7T7r=) ■ (3-43) 



2V47T, 

This CFT differs from that analyzed in the previous subsection for the exchange of the 
role of the two fundamental exponentials. 

According to the CFT defined by the bare lagrangian ( ft. 42] ) , the anomalous dimensions 
of the primary fields corresponding to the operators e Q< ^ are given by 

. , . a 2 aQ + (q) 

07T V47T 



and for the fundamental exponentials we have in this case 

A + (*-) =-!-£. <3 ' 45> 

Repeating the same kind of computations of the previous subsection, it is easy to check 
that these expressions coincide, at order 0(g 4 ), with the conformal dimensions extracted 
from the ultraviolet behaviour of the correlators < T zz (z, z) $±(0) >. 

For what concerns the central charge, since it depends quadratically on Q±, its value 
is given as before by the self-dual function ( ft.31| ). Analogous computations for the second 
moment of the G computed in terms of the FF of <&_(x,B) (which is the dual operator 
of < &+(x, £>)), lead therefore to the same results of Table 1 (see also Fig. 1). 



3.4.3 General case 

We are now able to give the conformal dimension A(a) of the exponential operator e aip 
and the central charge of the CFT reached in the ultraviolet regime for generic value of 
the parameter a. 

Since A (a) is the coefficient of the most singular term obtained in the UV-limit of the 
correlation function < T zz (z,z) e av ^ >, and the FF of T zz (z, z) depends linearly on the 
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parameter a, the conformal dimension is given by 



A(a) = oA_(a) + (l-o)A+(a) = (3.46) 
a2 + « (l-2a)f^+ 9 ^ 



7T V47T V # 2V4vr 



with A±(a) defined in eqs. ( |3.33| ) and ( |3.44D . The coefficient in front of the linear term 



in a in eq. (p.46Q identifies the background charge and therefore the central charge of the 



CFT reached in the ultraviolet limit is given by 

(l-2a) 2 

= 1 + 24^-^. (3.47) 

We have checked the validity of this result with the computation of the central charge in 
terms of the first FF of the operator ( |3.25| ) . The comparison between them is shown in 



Fig. 3, varying a at fixed B. 

As last example of possible choices of the stress-energy tensor, observe that for a = 1/2 
we have c = 1, independent of the coupling constant. The corresponding expression of G 
is given by || 



e(x) = u v ' [e 9tp + e-^J . (3.48) 

This operator is manifestly self-dual and Z 2 -even. With this choice of a, the anomalous 
dimensions of the fundamental exponentials coincide, at lowest order, with their gaussian 
values 

A(±g) = ~f + o{g±) . (3.49) 

87T 

The check of the central charge obtained in this case has been done in the two-particle 
approximation in ref. @ and is reported here in table 2. 



4 The Bullough-Dodd massive boson 

Different ultraviolet scaling regimes induced by different choice of the stress-energy tensor 
can be easily discussed for another integrable theory involving an interacting bosonic field, 
the so-called Bullough-Dodd (BD) model. 

4.1 Basic Properties 



The Bullough-Dodd (BD) model is defined by the equation of motion p2, ^3[ ^4] 

i 

3A 



Ut p = Ig>( e - 2 ^-e^) . (4.1) 
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At the quantum level, the integrability of the model leads to the elasticity and factorization 
of the scattering processes. The spectrum of the model consists of a massive particle state 
A created by the elementary field (p. This particle appears as bound state of itself in the 
scattering process 

AxA->A-^AxA. (4.2) 
The corresponding S-matrix is given by JTT" 



S(J3,B) = /|(/?)/|_ § . (4.3) 

Here 

fM = tanh ^ + mx) , (4.4) 
tanh — mx) ' 

and the coupling constant dependence of the model is encoded into the function 

H(A) = . (4.5) 

47T 

Like the Sinh-Gordon model, the S-matrix of the BD model is invariant under the mapping 

5(A) _> 2 -0(A) , (4.6) 

i.e. under the weak-strong coupling constant duality A — > 4tt/\. 

The minimal part of the S-matrix, i.e. the term contains the physical pole 

(3 = 27ri/3 of the bound state and, as matter of fact, it coincides with the S-matrix of 
the Yang-Lee model [23]. Taking into account the coupling constant dependence of the 



S-matrix, the residue at the pole is given by 



_ tan(^) tan (2 



6 

tan [if - 2p) tan(^f + f 



(4.7) 



This function, that corresponds to the three-particle vertex on mass-shell, vanishes for 
B = and B = 2 (which are the free theory limits) with the corresponding scattering 
amplitude S = 1. However, it also vanishes at the self-dual point B = 1, with the 
corresponding scattering amplitude S([3) = f-2/3- This coincides with the S-matrix of the 
Sinh-Gordon model computed at B = 2/3. As analyzed in [2B|, this equality between the 



S-matrices of the two models implies that at the self-dual point the BD model dynamically 
develops a ^-symmetry which is a non-perturbative property of the model. 



4.2 Form Factors 

Taking into account the bound state pole in the two-particle channel at /3y = 2m/ 3 and 
the one-particle pole in the three-particle channel at (3^ = in, the general form factors of 
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the BD model can be parameterized as 

i<j [ x i + Xj){uJXi + Xj){0J 1 X i + Xj) 

where we have introduced the variables 

Xi = e ft , uj = e in/3 . (4.9) 

FmWiiP) is an analytic function without zeros and poles in the physical sheet, whose explicit 
expression is given by 



jffiftB) = n 



fc=0 



r 


[k + l + i) 


r 


[k + l + i) 


r(k + l + i) 


r 


[k + l + i) 


r 


[k + l + i) 


r(k + l + i) 



(4.10) 



^(k+i-f+i)r{k+i+i+i)T(k + i-i + i)r( k+ i + i + i 



r(* + | + f + S)r(Jfe + i-f + S)r(* + i + f + ||)r(jfe + |-f + J|) 

{(3 = in — (3). Its normalization is fixed by requiring the asymptotic behaviour 



Notice that at the self-dual point 3 = 1, the above function coincides with the F^ n U3, |j 
of the Sinh- Gordon model, i.e. 

*ffi/M) = F ^(^l) ■ ( 4 - 12 ) 

The functions Q n (xi, ■ ■ ■ , x n ) are symmetric polynomials in the variables x±, . . . , x n . Using 
the functional relations satisfied by F min (f3, B) 

sinh (3 (sinh (3 + sinh ^ 

^S(^+Afi)i^£(/3,fi) ^ 



(sinh (3 + sinh ^) (sinh (3 + sinh ^ 



-B) 



-j -J cosh (3 + cosh — — - 

the kinematical and bound state residue conditions give rise to the following recursive 
equations satisfied by the functions Q n (xi, . . . , x n ) [p6fl 

1 

( -Q %2i ■ ■ ■ i %n) rpBD ( • ^(^"' "^l' "^2; ■ ■ ■ > x n)Qni. x li %2i • • • j x n) ) 

(4.14) 
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where 



U(x, Xl ,...,x n ) = 2 £ {-l)*^^*+-"^a%V£K..a$? (4.15) 

ki,...,ke=0 



x sm 

and 



| [2(fc 2 + A; 4 - fci - fc 3 ) + B(h + h-k A - fc 5 )] 



a/3 

Q n+2 (^£, u; _1 x, xi, . . . , x n ) = Bp t^- x r(#) x 3 L>(x, Xi, . . . , x n ) Q„+i(x, x : , . . . , x n ) , 

-^min ^ 3 ) 

(4.16) 

where 

n 

D(x,x 1 ,...,x n ) = H(x + Xi){xu 2+B + Xi)(xu~ B - 2 + Xi) (4.17) 



n 



E 



„3n-(fei+fc 2 +fc 6 ) , ,(2+S)(fe 2 -fc 3 ) An) (n) (n) 
x u °k x a k 2 a k :i ■ 



ki,k2,ke=0 

4.3 Cluster Operators and Fundamental Exponentials 

Unlike the Sinh-Gordon model, we do not know presently a close solution for the recursive 
equations satisfied by the FF of the BD model at generic value of the coupling constant. 
However, as it will become clear in the following, all we need for our consideration is the 
explicit computation of the first representative FF of the elementary field cp and of the 
so-called cluster operators V±(x). 

The elementary field (p is identified as that operator that creates one-particle state. 
Therefore its lowest matrix element is given by 

Ff =< 0\<p(0)\P >= ±= . (4.18) 

With such normalization, the next FF are given by 

F*(8) - - 3^ Sin ^ 2 + g ) U 19) 

and 

(tt F mil((3i ~ \ 2 a/2 

FM ' ft,ft) = [U (*i + xMeWxi + Xj )(e-^ Xl + Xj ) ) F^{in,Bf 3 * 

{2 Sin 2 I (2 + B) T 2 (B) (cOS I (2 + B) - l) (73(72(7! - (4.20) 

W{B) ((73(7? + af) + (sin 2 |(2 + B) T 2 (B) + W{B)) a 2 a 2 } , 
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where 

W(B) = 2V3 sin sin (^(2 - B)\ . (4.21) 

In terms of them, we can easily obtain the first FF of the operator Uip 

F °<p = _ m 2^zi F ^ . (4.22) 

In order to define the FF of the two fundamental exponential operators $1 (x, B) = 
e \(p{x) anc | <j) 2 ( X)J g) = e -2Ai^(a) Q £ foe Bullough-Dodd model, let us consider initially the 
definition of the cluster operators V±(x, B). As for the Sinh-Gordon model, we are looking 
for these operators in the class of FF which are asymptotically constant for x\ — > 00 and 
proportional to the invariant combination of symmetric polynomials 0\O n -\ for n > 2. 
The first representative of such FF are given in Appendix A. The important point is that 
all the higher FF obtained by solving the recursive equations will depend on the constants 
H\ and H 2 appearing in the equations ( A.l ) and ( |A.2j) , which play the role of arbitrary 



initial conditions of the recursive structure. Their relative value can be fixed though, if 
we require that the above FF satisfy an additional condition, i.e. the cluster property 



lim F k+l (ft + A, . . . , ft + A, ft +1 , . . . , (3 k+ i) = F k (ft, . . . , ft) Ft (ft +1 , . . . , f3 k+l ) 



(4.23) 



with Fq = 1. In this case we have 



H?(B) = . 1 \ - sin (^{B + 2)) T{B)± Jsin 2 (^(B + 2))T^B) + AW(B) 



Hf{B) = {Htf{B) . (4.24) 



With this choice of Hf(B) and Hf{B), the infinite tower of FF with JO|) and (|A~2 



as 



first representatives, define two cluster operators V±(x,B). By construction, the matrix 
elements of such operators are invariant under the duality transformation B — > 2 — B. 
Also in this case, we conjecture that the fundamental exponential operators are given by 

$x(x, B) = = 6(1 - B) V+{x, B) + 0(B - 1) V_(x, B) , 

$ 2 {x,B) = ee 6(1- B)V-(x,B) + 6(B-1) V+(x,B) . { ' ' 

This definition is in agreement with the perturbative analysis of the matrix elements of 
the two exponential operators e Xip ^ and e~ 2Xtp ( x \ Concerning their properties under the 
duality mapping, as far as B 7^ 1, the fundamental exponentials are mapped each into the 
other under the mapping B — > 2 — B, i.e. 

* lt 2(x,B) = <f> 2)1 (x,2-B) . (4.26) 

However, this mapping becomes degenerate at the self-dual point B = 1 where, similarly 
to the Sinh-Gordon model, the two operators $i(a;, B) and $2(2, B) collapse into a single 
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operator Moreover, as already noticed in [^J, at the self-dual point a Z 2 -symmetry 

is dynamically implemented in the BD model. Due to the fact that at B = 1 the three- 
particle vertex T(B) vanishes and Hi = —H{, the resulting operator will have 
non-zero matrix elements only on the 2n particle states. Its FF are entirely expressed in 
terms of the FF of the Sinh-Gordon model at B — 2/3 

^-(,(H)) 2 "Q„(1)I1^|^, (4.27) 

with M(B) and fi(B) defined in eqs. O and ([TTU]). 

Comparing the form factors of Oip and the form factors of the fundamental exponen- 
tials the quantum equation of motion can be cast in the form 



Uip(x) = F ^{ m ^ (6(1 -B) - 9(B - 1)) (e- 2X ^ - e x ^) 

2^2 \ sin 2 (f(£ + 2))r 2 (£)+4^(£) V 1 



2 

m 



F B P(m B) 

mm{ ' ' (V-(x,B) -V+(x,B)) . (4.28) 



2V2 \ sin 2 (f (B + 2)) T 2 (B) + 4W(B) 
4.4 Class of Stress-Energy Tensor 

The most general expression of the trace of the stress-energy tensor compatible with 
the (quantum) equation of motion of the BD model can be expressed in terms of the 
fundamental exponentials as 

Q(x) = F®(B)(a^ 1 (x,B) + (1 - a)$ 2 (x,B)) , (4.29) 

where F®(B) is its vacuum expectation value 

tf(B) = ^ , (4,0) 



(as computed by the Thermodynamical Bethe Ansatz p0|), whereas a is a free parameter. 
Varying the value of a, we may reach different ultraviolet limit of the BD model. Before 
considering the general case, let us analyze separately the two cases a = 1 and a = 0. 

4.4.1 The case a = 1 

For this value of a, the trace of the stress-energy tensor is given entirely by the operator 
$i(x,B) and therefore we expect that the ultraviolet behaviour will be described by a 
CFT with bare action given by 



d x 



2 W) 6A 2 e 



(4.31) 
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and background charge 



7T 



Q 2 (A) = -[ 2 y 



A 



(4.32) 



Using eq. 



the corresponding central charge is given by 
c(B) = 1 + 12Q' 2 2 (\) = 1 + 12 :1 ~ 5 6 



AB 



2-B 



(4.33) 



This is confirmed by the computation of the central charge in terms of the c-theorem by 
using the FF of the operator <$>i(x,B) which defines in this case the trace of the stress- 
energy tensor. The result of this computation is reported in Table 3 (see also Fig. 2) and 
the sum rule turns out to be saturated with high percentage of precision by using just the 
first two FF of $i(x,B). 

According to the CFT ( 4.31f) , the conformal dimensions of the primary fields e aip are 
given by 

«Q 2 (A) 



A 2 (a) = h 

1 ; 8tt V47T 

and for the fundamental exponential operators of the BD model we have 



(4.34) 



A 2 ($ 2 ) 
A 2 ($i) 



— A 2 

8vr 



(4.35) 



It is easy to see that these expressions are in agreement with those extracted by looking at 
the short-distance behaviour of the correlators < T zz (z, z) 3>i(0) > and < T zz (z, z) $ 2 (0) >. 
The computation are similar to that of the Sinh-Gordon model, the only difference being 
the perturbative expansion 



exp 



dx sinh (jf \ sinh (|(2 - Bn sinh | cosh ■ 



io x 
'1 1 

7T 



sinh x 



e7S) T + o(A4 » 



(4.36) 



and therefore we will not repeat them here. 



4.4.2 The case a = 

Since the trace of the stress-energy tensor is given in this case by the operator $ 2 (a;, B), 
the ultraviolet limit will be ruled by a CFT with a bare action given by 

,2 



Si 



d 2 x 



m 
3A 1 



(4.37) 
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and a background charge given by 

Qi(A) 



+ 



A 



V A 2V47T, 
For the corresponding value of the central charge we have 



c(B) 



1+ 12Q 2 (A) 

,'2-B IB 
1 + 12 — — + 



B 42-S 

whereas for the conformal dimension of the primary operators e aip 

Ax(a) = + .. 
Hence, for the fundamental exponential operators we have in this case 



Ai($0 
Ai($ 2 ) 



-2 - —A 2 

47T 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



These conformal data are again confirmed by the FF approach, as shown for instance in 
Fig. 2 for the central charge. 



4.4.3 The general case 

It is now easy to write down the conformal dimension A (a) of the exponential operator 
e atp and the central charge of the CFT reached in the ultraviolet regime for generic value 
of the parameter a appearing in the definition of the stress-energy tensor of the BD model. 
The argument is similar to that already employed in the Sinh-Gordon model. Since A(a) 
is the coefficient of the most singular term obtained in the UV-limit of the correlation 
function < T zz (z, z) e a ^ 0) > and the FF of T zz (z, z) depends linearly on the parameter a, 
the conformal dimension is given by 



A(a) = aA 2 (a) + (1 

a 2 a 
8n J^n 



a) Ai(a) 
\l-B 



(4.42) 



B 



3 

- a 
2 



+ 



B 



B 



3 

- a 
2 



where Ai(a) and A 2 (a) are given in eq. Q4.34D and Q4.40D respectively. The linear term in 
a in ( |4.42| ) identifies the background charge of the corresponding Coulomb gas. Therefore 
the central charge of the CFT reached in the ultraviolet limit is given by 



c = 1 + 3 



(B + Qa-Af 
B{2 - B) 



(4.43) 
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The check of this formula in terms of the FF approach is shown in Fig. 3. 
Observe that, with the choice 

4 - B 

a = , (4.44) 

we have identically c = 1, a result that is confirmed by the FF approach within the usual 
accuracy of few percents. The corresponding trace of the stress-energy tensor is given by 

G(x) = Ff(B) fi^e A ^) + m e -2M*)\ . (4.45) 
V 6 6 / 

This operator is manifestly self-dual. In the limit A — > 0, it reduces to 

e(ar) = ^( 2 e A ^) + e- 2A ^)) , (4.46) 
which is the classical expression of Q(x) for the Bullough-Dodd model. 



5 Conclusions 

The ultraviolet behaviour of a two-dimensional QFT is generally characterized by a scaling 
behaviour described by a CFT. The main features of a CFT are encoded in the definition 
of the stress-energy tensor T^ix). As shown by the form factor approach, associated 
to an on-shell dynamics, there is a one-parameter family of possible operators T /lu (x) 
that induces different scaling behaviour of the massive theory in the ultraviolet limit. 
In light of their simple spectrum, we have analyzed in detail the Sinh-Gordon and the 
Bullough-Dodd Models, computing the relevant CFT data (central charge and conformal 
dimensions) in terms of the FF of the fundamental exponential operators. 

The ultraviolet properties of these models are strictly related to the duality symmetry 
of their S- matrix. This symmetry has far-reaching consequences. In fact, since the FF 
are computed in terms of the 5-matrix, this symmetry may also be extended off-shell. In 
particular, it gives rise to a bidimensional representation in the space of the fundamental 
exponential operators of both theories. Moreover, from the self-duality of the two theories, 
we have an identification of the two operators at the self-dual point. This remarkable 
property is based on our definitions (|3.20|) and ([4.25 ) of the fundamental exponentials 



of the two theories in terms of their respective cluster operators. The validity of these 
formulas has been checked by analyzing the ultraviolet behaviour of the massive theories. 
An important difference between the BD and the SG models is in the expressions of the 
central charge associated to their fundamental exponentials. In fact, whereas in the SG 
model we have the self-dual function ( 3.3 1| ) which is common to both the exponentials 



&±(x), in the BD model we have the different expressions (|4. 33|) and (|4.39|) , related, 
however by duality 

c(a= 0,B) = c(a = 1,2- B) . (5.1) 
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It would be interesting to extend our results to other Afnne Toda Field Theories and to 
study in more detail the interplay between the massive and conformal data of the models. 
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Appendix A 



Solving the recursive equations ( |4.14| ) and fl4.16|) , the first FF which have the properties to 



be asymptotically constant and proportional to the invariant combination o\o n -\ (n > 2) 
are given by 

and 

F3{Pi,(h,Pa) ( x . + Xj )( UXi + Xj )(u- 1 x i +Xj) ) F^P(iir,B) e7l ' T2 * 

{qx(B)al + q 2 (B)a 3 a 2 a l - (A.2) 
H X W{B) (a 3 af + al) + q*{B)ala{} , 



where 



Ql (B) = 2r(B)sin£(2 + B) 

o 

q 2 (B) = T(B)sm^(2 + B) 
o 



cos^(2 + £) -1 



vr _ ~ H 2 



r(B)H lS m-(2 + B) +y/F*&(in)-£ 



Y{B)H X sin |(2 + B) - y/F*P(in)H 2 (cos |(2 + B) - - 
(6(B) = ff^(B) - y/F^(in)^T(B) sin^(2 + B) . (A.3) 
.Hi and i?2 are two arbitrary parameters. 
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Table Caption 



Table 1 . The first two-particle term entering the sum rule of the c-theorem for the Sinh- 
Gordon model with the choice a = 1 (second column) compared with the central 
charge ( ggg) of CFT. 



Table 2 . The first two-particle term entering the sum rule of the c-theorem for the 
Sinh-Gordon model with the choice a = 1/2. It must be compared with the claimed 
free boson Cuv — 1 UV-behaviour. 

Table 3 . The first two-particle term entering the sum rule of the c-theorem for the 
Bullough-Dodd model with the choice a = 1 (second column) compared with the 
central charge ( |4.33| ) of CFT. 



28 



B 




/^Liouv 


1 

10 


127.28994 


127.31579 


1 

5 


67.61695 


67.66667 


3 
10 


47.98763 


48.05882 


2 
5 


38.40998 


38.5 


1 

2 


32.89395 


33. 


3 
5 


29.45222 


29.57143 


7 

10 


27.24418 


27.37363 


4 

5 


25.86323 


26. 


9 
10 


25.10126 


25.24242 


1 


24.85738 


25. 



Table 1 



29 



B 




1 

500 


0.9999995 


1 
100 


0.9999878 


1 

10 


0.9989538 


3 
10 


0.9931954 


2 
5 


0.9897087 


1 

2 


0.9863354 


2 
3 


0.9815944 


7 

10 


0.9808312 


4 

5 


0.9789824 


1 


0.9774634 



Table 2 



30 



B 




/^Liouv 


1 

10 


70.63001 


70.63158 


1 

5 


41.32883 


41.33333 


3 
10 


32.10844 


32.11765 


2 
5 


27.98391 


28. 


1 

2 


25.97441 


26. 


3 
5 


25.10474 


25.14286 


7 

10 


24.97886 


25.03297 


4 

5 


25.42607 


25.5 


9 
10 


26.38691 


26.48485 


1 


27.87364 


28. 


11 
10 


29.96195 


30.12121 


6 
5 


32.80360 


33. 


13 
10 


36.66406 


36.90110 


7 
5 


42.00619 


42.28571 


3 
2 


49.67928 


50. 


8 
5 


61.39520 


61.75 


17 
10 


81.15833 


81.52941 


9 
5 


120.98370 


121.33333 


19 
10 


240.90584 


241.15789 



Table 3 



31 



Figure Caption 



Figure 1 . The first two-particle term entering the sum rule of the c-theorem for the 
Sinh-Gordon model with the choice a = 1,0 (dots), compared with the selfdual 
central charge ( |3.31| ) of CFT (solid line). 

Figure 2 . The first two-particle term entering the sum rule of the c-theorem for the 
Bullough-Dodd model with the choice a = 1,0 (dots and crosses, resp.), compared 
with the non selfdual central charges fl4.33|) , ( |4.39|) of CFT (solid thin line and solid 
thick line, resp.). 

Figure 3 . The first two-particle term entering the sum rule of the c-theorem for the 
Sinh-Gordon model and the Bullough-Dodd model with the coupling constant fixed 
at B — 1/2 for different values of a (crosses and dots, resp.), compared with the 
central charges ( |3.47| ), ( |4.43| ) of CFT (solid thin line and solid thick line, resp.). 
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